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Summary. — Electron-positron plasma is believed to play imporant role both 
in the early Universe and in sources of Gamma- Ray Bursts (GRBs). We focus on 
analogy and difference between physical conditions of electron-positron plasma in 
the early Universe and in sources of GRBs. We discuss a) dynamical differences, 
namely thermal acceleration of the outflow in GRB sources vs cosmological deceler- 
ation; b) nuclear composition differences as synthesis of light elements in the early 
Universe and possible destruction of heavy elements in GRB plasma; c) different 
physical conditions during last scattering of photons by electrons. Only during the 
acceleration phase of the optically thick electron-positron plasma comoving observer 
may find it similar to the early Universe. This similarity breaks down during the 
coasting phase. Reprocessing of nuclear abundances may likely take place in GRB 
sources. Heavy nuclear elements are then destroyed, resulting mainly in protons 
with small admixture of helium. Unlike the primordial plasma which recombines 
to form neutral hydrogen, and emits the Cosmic Microwave Background Radiation, 
GRB plasma does not cool down enough to recombine. 
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1. Introduction 

Electron-positron plasmas are discussed in connection with astrophysical phenomena 
such as Galactic Center, microquasars, Gamma- Ray Bursts (GRBs), as well as labora- 
tory experiments with high power lasers, for details see [1]. According to the standard 
cosmological model, such plasma existed also in the early Universe. It is naturally char- 
acterized by the energy scale given by the electron rest mass energy, 511 kcV. It is 
interesting that at the epoch when Universe had this temperature, several important 
phenomena took place almost contemporarily: electron-positron pair annihilation, the 
Big Bang Nucleosynthesis (BBN) and neutrino decoupling. 

Electron-positron plasma also is thought to play an essential role in GRB sources, 
where simple estimates for the initial temperature give values in MeV region. Such 
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plasma is energy dominated and optically thick due to both Compton scattering and 
electron-positron pair creation, and relaxes to thermal equilibrium on a time scale less 
than 10 -11 sec, see [2]. The latter condition results in self-accelerated expansion of the 
plasma until it becomes either transparent or matter dominated. 

In the literature there have been several qualitative arguments mentioning possible 
similarities between electron-positron plasmas in the early Universe and in GRB sources. 
However, until now there is no dedicated study which draws analogies and differences be- 
tween these two cases. This paper aims in confronting dynamics and physical conditions 
in both cases. 

2. General equations 

The framework which describes electron-positron plasma both in cosmology and in 
GRB sources is General Relativity. Both dynamics of expansion of the Universe, and 
the process of energy release in the source of GRB should be considered within that 
framework. Hydrodynamic expansion of GRB sources may, however, be studied within 
much simplier formalism of Special Relativity. 

We start with Einstein equations 

1 8ttG 

(1) R^v — ~^3llvR — ^ 4 Tpv , 

where R^, g^ v and T^ v are respectively Ricci, metric and energy-momentum tensors, 
G is Newton's constant, c is the speed of light, and the energy-momentum conservation, 
following from (1) 

(2) (T, v ), v = 4= diV 7l" - T ^ = 0> 



dx v 

where are Cristoffcl symbols and g is determinant of the metric tensor. We assume 
for the energy-momentum tensor 

(3) T» v = pg" v +ljU^U u , 

where J7 M , is four- velocity, to = p + p is proper enthalpy, p is proper pressure and p is 
proper energy density. 

When plasma is optically thick, radiation is trapped in it and entropy conservation 
applies. It may be obtained multiplying (2) by four-velocity 

(4) —U^ (Tju ") .„ = U^p-fj, + uiU^-fj, = 0. 
Using the second law of thermodynamics 

(5) d(^)=T d r-) + 1 -dp, 

where a = uj/T is proper entropy density, T is temperature, one may rewrite (4) as 

(6) {oU»), ll = U»a. ll + aU». lx = Q. 
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Baryon number conservation equation has exactly the same form 



(7) 



(n[/'% = l^7i ;M + n[/% = 0. 



Now recalling that U 11 -^ = g~ and U^ ;tl = dlnV/dr, where V is comoving volume, 
r is the proper time, from (4) and (7) we get 



(8) 



dp + udlnV = 0, dlnn + dlnV = 0, 



Finally, introducing the thermal index 7=1 + - restricted by the inequality 1 < 7 < 4/3 
we obtain from (4) the following scaling laws 



(9) 



pV 1 = const, nV = const. 



Both these conservations laws are valid for the early Universe and GRB plasmas. 

One can obtain the corresponding scaling laws for comoving temperature by splitting 
the total energy density into nonrelativistic (with 7 = 1) and ultrarclativistic (with 

e, where m is the mass of particlcs( 1 ), e is proper 



7 = 4/3) parts with p — > nmc 2 
internal ener; 
and (6) gives 



internal energy density. The entropy of the ultrarelativistic component is then a = | A 



(10) 



eV 
~T 



const. 



For e ^> nmc? , which is the energy dominance condition, internal energy plays dynamical 
role by influencing the laws of expansion. For e -C nmc 2 , which is the matter dominance 
condition, internal energy does not play any dynamical role, but determines the scaling 
law of the temperature. In order to understand the dynamics of thermodynamic quan- 
tities in both early Universe and in GRBs, one should write down the corresponding 
equations of motion. 

2T. Early Universe. - For the description of the early Universe we take the Robertson- 
Walker metric with the interval 



(11) 



ds 2 = -c 2 dt 2 + a 2 (i) 



dr 2 



1 - fcr 



+ r 2 d-d 2 + r 2 sin 2 ddip 1 



where a (t) is the scale factor and k = 0,±1 stands for the spatial curvature. In ho- 
mogeneous and isotropic space described by (11), Einstein equations (1) are reduced to 



( x ) Nonrelativistic component is represented by baryons. For simplicity we assume only one 
sort of baryons, say protons, having mass m. Ultrarelativistic component is represented by 
photons and electron-positron pairs. 
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Firedmann equations together with the continuity equation 
,„ f da\ 2 , , 8nG 9 

< 12 > UJ + c2fc =^ 2 - 

. . „ d 2 a ( da\ 2 2 , 87rG 2 

(13) 2a ^+Uj +c 2 ^-— P a 2 , 

dp 3 da , 

where a is the scale factor. Notice, that only two equations in the system above are 
independent. The continuity equation (14) follows from the Einstein equations (12) and 
(13) as the energy conservation. In fact, (14) may be also obtained from the entropy 
conservation (4). The comoving volume in Friedmann's Universe scales with a as V = a 3 , 
so (14) and the first equality in (9) are equivalent. 

On the radiation dominated stage of the Universe expansion one has 

(15) p cx U~ 4/3 oc a" 4 , iKxF'aa- 3 , 
while on the matter dominated stage 

(16) p oc n oc V^ 1 oc a~ 3 . 

Entropy conservation (10) leads to the unique temperature dependence on the scale factor 

(17) TocV~ 1/3 



oc a 1 . 



The corresponding time dependence of thermodynamic quantities may be obtained 
from solutions of Fricdmann equation (12) and continuity equation (14), see e.g. [3]. 

2'2. GRBs. - Different situation takes place for the sources of GRBs. Assuming 
spherical symmetry for the case of GRB the interval( 2 ) is 

(18) ds 2 = -c 2 dt 2 + dr 2 + r 2 dd 2 + r 2 sin 2 Mtp 2 . 

Optically thick to Compton scattering and pair production electron-positron plasma in 
GRB sources is radiation dominated. Its equations of motion follow from the energy- 
momentum conservation law (2) and baryon number conservation law (7). Initially 
plasma expands with acceleration driven by the radiative pressure. 

In spherically symmetric case the number conservation equation (7) is 



( 2 ) General Relativity effects may be included by taking Schwarzschild or Kerr-Newman met- 
ric. However, we are interested in optically thick plasma which expands with acceleration and 
propagates far from its source, where the spatial curvature effects may be neglected. For this 
reason we simplify the treatment and adopt a spatially flat metric. 
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Integrating this equation over the volume from certain (t) to r e (t) which we assume to 
be comoving with the fluid dri Jp = /3(ri,t), = j3(r e ,t), and ignoring a factor 4n we 

have 



(20) J ^^r 2 dr + (r 2 n^f^~l) dr = 

Q f df dv ' 

— J (nT) r 2 dr - ~^n(r e , t)T(r e ,t)r 2 e + -£n(n, t)T(n, t)rf + 

+r 2 e 7i(r e ,t)^{r e ,t) - 1 - r 2 n{ n , t) ^/T 2 { n , t) - 1 = 

T e 

d f 

= — / (nT) r 2 dr = 0, 

Since we deal with arbitrary comoving boundaries, this means that the total number of 
particles integrated over all differential shells is conserved 



j 1,1 1 

(21) 7V = 4tt J 



R(t) 

r / nTr 2 dr 
" 



where R(t) is the external radius of the shell. 

Following [4] one can transform (19) from the variables (t, r) to the new variables 
(s = t — r, r) and then show that 

For ultrarelativistic expansion velocity r>l, the RHS in (22) tends to zero, and then 
the number of particles in each differential shell between the boundaries 7"j(i) and r e (i) 
is also conserved with a good approximation, i.e. 



(23) dN = &imTr 2 dr « const. 

Relations (21) and (23) then imply 



(24) 4tt / (nFr 2 ) dr = 4tt [n(r, t)T(r, t)r 2 ] J dr 

Ti Ti 

= An (nTr 2 ) A ps const, 

where the first argument of functions n(r, t) and T(r, t) is restricted to the interval < 
r < r e , and consequently A = r e — ri ~ const. Taking into account that r,(i) and r e (f) 
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are arbitrary, this means that ultrarelativistically expanding shell preserves its width 
measured in the laboratory reference frame. This fact has been used in [5] and referred 
there as the constant thickness approximation. 

The volume element measured in the laboratory reference frame is dV = 4irr 2 dr, 
while the volume element measured in the reference frame comoving with the shell is 
dV = 4irTr 2 dr. Comoving volume of the expanding ultrarclativistic shell with T ~ const 
will be 

(25) V = 4tiT / r 2 dr ~ 4irTr 2 A. 

Jr-A 

Then we rewrite the conservation equations (9) as 

(26) p~Tr 2 = const, nTr 2 = const, 

Unlike the early Universe, where both energy and entropy conservations reduce to (14), 
in the case of GRBs the energy conservation is a separate equation coming from the 
zeroth component of (2) as 

(27) (T v ) , v = ljU U v - v + U v (uU ). u = 0. 
which independently of 7 gives 

(28) pY 2 r 2 = const. 
From (26) and (28) we then find 

2(7-1) 2_ 27 

(29) r oc r 2 -~t , n oc r 2 -y , p oc r 2 -t . 

For the ultrarelativistic equation of state with 7 = 4/3 we immediately obtain 

(30) r a r, n oc r~ 3 , poc?*~ 4 . 

Taking into account that the relation between the comoving and the physical coordinates 
in cosmology is given by the scale factor a, it follows from (30) that both energy density 
and baryonic number density behave as in the radiation dominated Universe, see (15). 
This analogy between the GRB source and the Friedmann Universe is noticed by [6] , [4] . 

In the presence of baryons as the pressure decreases, plasma becomes matter domi- 
nated and expansion velocity saturates. Hence for the nonrelativistic equation of state 
with 7=1 different scaling laws come out 

(31) r = const, n oc r~ 2 , pocr" 2 . 

Transition between the two regimes (30) and (31) occurs at the radius R c = i? _1 i? , 
where i?o is initial size of plasma. 

Therefore, one may reach the conclusion that for comoving observer the radiation- 
dominated plasma looks indistinguishable from a portion of radiation-dominated Uni- 
verse. However, this is true only in the absence of pressure gradients. Strong gradients 
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are likely present in GRB sources, and they should produce local acceleration in the 
radiation-dominated electron-positron plasma, making it distinct from the early Uni- 
verse, where matter inhomogeneitics arc known to be weak. 

It is easy to get from (26) and (28) for internal energy density and temperature 



(32) eocr~ 4 , TocT 1 , R < r < R c , 
and 

(33) eocr" 8 / 3 , T«r 2/3 , R c < r < R, 



where Rt r is the radius at which the outflow becomes transparent. The outflow may 
become transparent for photons also at the acceleration phase, provided that Rt r < R c - 
For instance, a pure electron-positron plasma gets transparent at the acceleration phase. 

3. Heavy elements 

Cosmological nucleosynthesis is a well established branch of cosmology. Classical 
computations made in the middle of the XXth century revealed that heavy elements 
cannot be built in the early Universe. Hydrogen and helium contribute approximately 
3/4 and 1/4, leaving some room, much less than 1 per cent for deuterium, tritium and 
lithium. All the heavier elements must have been produced in stars. 

Some of these stars, as indicated by observations, end their life as progenitors of 
GRBs. For this reason it is likely that initially in the source of GRBs elements heavier 
than hydrogen are present. In this section we consider chemical evolution of plasma in 
the sources of GRBs. 

Assume that in the source of a GRB the amount of energy E$ is released in the volume 
with linear size Rq during the time At, making this region optically thick to Compton 
scattering and pair production. The amount of baryons which may be present as well is 
parametrized by 



2 



(34) B 



Mc 

Mc 2 
L 



, At < Rq/c, 
, At > Rq/c, 



where L = dE/dt is the luminosity, M = dM/dt is the mass ejection rate and M is 
total baryonic mass. Ultrarclativistic outflow is generated through thermal acceleration 
of baryons by the radiative pressure if plasma is initially energy dominated, i.e. 

(35) 



In the case of instant energy release with time interval At ~ Ro/c initial temperature 
in the source of GRB may be estimated neglecting the baryonic contribution, provided 
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(35) is satisfied as 



( 



3£o 



) 



1/4 



(36) 



To 



6.5£;^ 4 i?8 3/4 MeV, 



where a = Agsb/c, usb is the Stefan-Boltzmann constant and the last value is obtained 
by substituting numerical values for Eq = 10 54 i?54 erg and Rq = 10 8 i?s cm. 

As it has been shown in [7] for temperatures above 1 MeV even low density plasma 
with density n = 10 18 cm~ 3 quite quickly destroys all heavier nuclei, and the final state 
contains just protons and neutrons and some small traces of Deuterium and 4 He. The 
timescale of this process (~ 10 -2 sec for Tq = 1 MeV) strongly depends on temperature, 
but the rates of almost all reactions increase with temperature, and correspondingly the 
abundances of nuclei evolve much faster. Therefore, nuclei disintegration is fast enough 
to occur before plasma starts to expand and cool on the timescale Rq / c. 

During early stages of plasma expansion its temperature decreases in the same way as 
it happened in the early Universe. Therefore similar synthesis of light elements to BBN 
occurs also in sources of GRBs. Most important is, however, another similarity with the 
early Universe: it is well known that practically all free neutrons have been captured 
into elements heavier than hydrogen. So we do not expect dynamically important free 
neutrons present in GRB plasma after it started to expand and cool down unless they 
are engulfed by the expanding plasma later. The role of such free neutrons have been 
considered in the literature, see e.g. [8] and [9]. 

4. Recombination 

On the radiation dominated phase both in the early Universe and in GRB plasma 
entropy conservation (4) results in decrease of temperature. When the comoving tem- 
perature decreases below the hydrogen ionization energy, Ei = 13.6 eV, the formation of 
neutral hydrogen is expected. 

4T. Early Universe. - In the early Universe, after the BBN epoch and electron- 
positron annihilation, cosmological plasma consists of fully ionized hydrogen, helium and 
small admixture of other light elements. The temperature continues to decrease until 
it gets sufficiently low to allow formation of neutral atoms: that is the moment in the 
cosmic history where the formation of the Cosmic Microwave Background Radiation 
(CMB) happens. 

The theory of cosmological recombination of hydrogen, based on three level approxi- 
mation, has been developed in [10] and [11] in the late 60s. The only modification that 
such theory undergone in the later years is the account for dark matter and addition of 
more levels to the model, currently about 300. There is a basic difference with respect 
to the equilibrium recombination essentially by the process e + p <-> H + 7, described by 
the Saha equation 



where gi are statistical weights, h is Planck's constant. This difference is due to the 
presence of the 2p quantum level, which produces Ly-a photons. The absorption of 
such photons is very strong. However, ionization from the 2p level requires only Bj/4. 



(37) 
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Therefore the formation of neutral hydrogen proceeds through the 2s — Is transition in 
the presence of abundant Ly-a photons. 

In fact, the early Universe would become transparent for radiation even if formation 
of hydrogen would have been forbidden, see e.g. [12]. The optical depth to Thomson 
scattering is 



*0 

aTTibcdt 



(38) 4xl0" 2 ^p^ n A + fi™(l + ^ 

" "177. 



1/2 



where <jt is the Thomson cross section, fli = Pi/p c , Pc = 3H 2 c 2 /8irG, H = lOOh km 
s/Mpc and 6, m, A stand for, respectively baryons, dark matter and cosmological constant 
contributions to the total energy density of the Universe. For large z we have 

r = 1, z» ~ 8.4fi- 2/3 0^ 3 /r 2/3 . 

For typical values il^h 2 ~ 0.02, TO ~ 0.3, and h ~ 0.7 we have ~ 60. At such redshift 
the Universe would be expected to become transparent to Thomson scattering. That is 
exactly what happens in plasma in GRB sources. Below we show that, unlike radiation- 
dominated cosmological expansion where comoving quantities also fulfill relations (32), 
the comoving temperature in GRB outflows remains always high enough to prevent 
recombination of hydrogen. 

4'2. GRBs. - During both acceleration and coasting phases the comoving tempera- 
ture decreases with radius, see (32) and (33). The optical depth to Compton scattering 
may be computed and the corresponding photospheric radius may be obtained, see [13] 
where ultrarelativistic outflows were analysed in details. In particular, comoving tem- 
perature at the photosphere decreases with the baryonic loading B both at acceleration 
and coasting phases. However, when the outflow reaches the radius R s = B~ 2 Ro, the 
comoving temperature becomes independent from B. In that regime the expression for 
the photospheric radius is 



(39) R tr 



( aE B N 

\47T77lpC 2 



where m p and a are proton mass and Thompson cross section, respectively. Indeed, using 
(36), (39), (32) and (33) in the case of instant energy release we have 



(40) T min = BT, 



R r N 2 / 3 



o , R 

\ 1/4 / x -1/3 



^ =Uw M {EoRo) 

Notice how extremely insensitive this value is with respect to the remaining parameters 
E and RqI Expressed in units of typical energy and size 

(42) T^ ) n ~42(i? 5 4 J R 8 )- 1/12 eV. 
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In the case of gradual energy release with At 3> Ro I c an d constant luminosity and mass 
ejection rate the initial temperature is 

/ L N 1/4 

(43) To " 



\ 16na S BRo 
and similar expression to (40) may be derived 

which may be rewritten, introducing L = 10 50 L5o erg/s and At = lAti s, as 

(45) T^cl7L- 1/12 Rl /6 At; 1/3 eV. 

Even if (45) appears to be less stringent that (42), they are both quite insensitive to 
initial parameters. As a result, even if the comoving temperature decreases very much 
compared to its initial value, typically on the order of MeV, at the photospheric radius it is 
always well above the characteristic temperature 0.3 cV at which recombination happens 
[14], thus preventing formation of neutral hydrogen. In fact, if such hydrogen would 
be formed the cross section of interaction of expanding particles with the circumburst 
medium would drastically decrease. As a consequence no afcrglow would be observed. 

A simplified way to look at this lower bound on the comoving temperature at the 
photosphere is to say that if a fraction e of solar mass is released in the volume hav- 
ing radius 8 solar Schwarzschild radii, then its minimum comoving temperature before 
transparency is 

(46) T^ n ^ 66 (eSy 1/12 eV, 
in the case of instant energy release and 

(47) T« ~ 2.8e- 1 /i2 (5 i/6 A i-i/4 eV> 

in the case of gradual energy release during time At±. Clearly in both cases 8 > 1, and 
likely e < 1. Notice, that while in the case of instant energy release the lower bound 
on temperature decreases with increasing 8, it instead increases in the case of gradual 
energy release. 

The baryon to photon ratio in GRB plasma like in cosmology is large. This ratio may 
be estimated as 

(48) ^ =^ c 1.8 x lO 5 ^ < £ }-\ 
n B m e B (e) 

where £?_ 2 = 10~ 2 B, and (e) is average photon energy in the source of GRB in units of 
electron rest mass energy. Thus the optical depth of electrons is much larger than the 
one of photons and it is given in [15] 

(49) !k = i ogA +^~!bL ) 

t 7 n e n B 
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where A is the Coulomb logarithm. It means that electrons are kept in equilibrium with 
photons when the latter already decoupled from them [16]. In other words, electrons 
are forced to keep the local temperature of photons. This may lead to efficient Comp- 
tonization of the photon flow when it is decoupled from plasma and is passing through 
electrons having locally different temperature. 

As soon as plasma gets collisionlcss, laboratory spectrum of photons, baryons and 
electrons is maintained. If it was thermal at decoupling it will remain so. This shows 
another difference with respect to cosmology, where energy of all particles decoupled 
from the thermal bath decreases due to the cosmological expansion. For that reason 
in cosmology only the shape of the spectrum is conserved with expansion, but not the 
temperature. 

Therefore, we have reached the conclusion that hydrogen recombination which is 
responsible for transparency of cosmological plasma does not occur in GRB plasma. This 
difference in physical conditions may result in deviations from black body spectrum, as 
observed in GRBs. Available studies of photospheric emission in GRBs in the literature 
show that deviations from the perfectly thermal spectrum come mostly from three effects: 
a) dynamical and ultrarelativistic character of plasma outflows and geometric effects 
[13]; b) "fuzzy photosphere" effect [17, 18] and c) possible dissipation mechanisms at the 
photosphere [9, 19, 20]. 

Recently we presented a theory of photospheric emission from relativistic outflows, 
see [13]. Assuming that the spectrum of radiation in the comoving reference frame is 
the perfect black body one, we have shown that the spectrum seen by a distant observer 
may be essentially nonthermal due to both geometric and dynamical special relativistic 
effects. The possibility that the spectrum of photospheric emission is nonthermal also in 
the comoving frame is under investigation. 

5. Conclusions 

Regarding the dynamical aspects, there is an apparent similarity between the electron- 
positron plasma in the early Universe and the one in GRB sources. For an observer 
comoving with the radiation-dominated plasma in GRB source it may look indistin- 
guishable from a portion of radiation-dominated Universe. However, this is true only in 
the absence of pressure gradients. Strong gradients are likely present in GRB sources, 
and they should produce local acceleration in the radiation-dominated electron-positron 
plasma, making it distinct from the early Universe, where matter inhomogeneities are 
known to be weak. 

There is also an apparent similarity with respect to the nucleosynthesis phenomenon. 
Given that the temperature reached in GRB sources, see Eq. (36), may be as high 
as several MeV, nuclear reactions are expected to operate on timcscalcs of 10~ 2 sec or 
shorter. That is on the order of magnitude of dynamical timescale of the GRB sources. 
It means that reprocessing of nuclear abundances may likely take place in GRB sources. 
Since observations imply that GRBs may originate from compact stellar objects elements 
heavier than helium are likely to be present in GRB sources. Such heavy elements are 
then destroyed, resulting mainly in protons with small admixture of helium. Thus, 
similarly to the early Universe, we do not expect dynamically important free neutrons 
present in GRB plasma after it started to expand and cool down unless they are engulfed 
by the expanding plasma later. 

Finally, there is an important difference between the electron-positron plasma in the 
early Universe and the one in GRB sources. We show in this paper that unlike the 
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primordial plasma which recombines to form neutral hydrogen, and emits the Cosmic 
Microwave Background Radiation, GRB plasma does not cool down enough to recombine. 
Therefore GRB plasma always becomes transparent due to Compton scattering. 
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